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Direct reconstruction of spherical harmonics from 
interferometer observations of the CMB polarization 
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ABSTRACT 

Interferometric observation of the CMB polarization can be expressed as a linear sum 
of spherical harmonic coefficients a±2,i m of the CMB polarization. The linear weight for 
a ±2,Zm depends on the observational configuration such as antenna pointing, baseline 
orientation, and spherical harmonic number l,m. Since an interferometer is sensitive 
over a finite range of multipoles, a±2,im m the range can be determined by fitting 
a±2,im for visibilities of various observational configurations. The formalism presented 
in this paper enables the determination of a±2,i m directly from spherical harmonic 
spaces without spherical harmonic transformation of pixcllized maps. The result of its 
application to a simulated observation is presented with the formalism. 
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1 INTRODUCTION 

The Cosmic Microwave Background (CMB) is expected to 
be linearly polarized by Thomson scattering at the last 
scattering surface and after re-ionization. T he CMB po- 
lariza tion has b een measured by the DA SI l|Leitch et alj 
120051), the CBI (iReadhead et alj[200i ) the BOOMERanG 
l|Montrov et al.ll2006t ), the CAPMAP dBarkats et alj|2005h 
and the WMAP satellite |Page et al.l 120061 ). Acharac- 
teristic signature imprinted on CMB polarization pro- 
vides valuable cosmological and astrophysical information. 
If the CMB anisotropy follows Gaussian distribution, the 
complete description of the CMB anisotr opy is provided 
through the angular po wer spectrum C; (|Dodelsonl 120031 ; 
iHinshaw and et alj|2007h . 

Interferometers offer more control of systematic ef- 
fects than traditional imaging systems, and have less E 
and B mode mixing (|Park and Nsd 12004 ) ■ With desirable 
features of interferometers, the CMB polarization mea- 
surement with interferometers is on- goi ng and planned in 
the experiments s u ch a s the DASI |Kovac et al.l 120021 : 
iLeitch et al.l |2002|. [20051). the Co smic Background Im- 
ager (CBI) (|Readhead et all 12004 ) and the Millimeter- 
wave Bolometric Inte rferometer (MBI) (|Tucker et al.|[20o3 : 
iKorotkov et al.ll2006l ). The usual procedure for the CMB 
analysis on the interferometer observation is to proceed to a 
statistical analysis such as maximum likelihood estimation 
of power spectra. In the power spectrum estimation by max- 
imum likelihood method, 0(Af 3 ) process should be carried 
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out repeatedly for iterative search. Since it becomes compu- 
tationally prohibitive wit h very large num ber of data, an un- 
biased hybrid estimator (Efstathiou 2006) proposes pseudo- 
Ci estimates at high multipoles, which requires the estima- 
tion of individual spherical harmonic coefficients. Spherical 
harmonic coefficients can be estimated also from mosaiced 
sky patches, which are reconstructed from interferometer 
observations via aperture synthesis. Due to flat sky approx- 
imation for each sky patch, the mosaiced sky map has dis- 
continuity on junctures of sky patches. For these reasons, 
we have investigated reconstructing the spherical harmonic 
coefficients of the CMB polarization directly from interfer- 
ometer observations in the complete context of a spherical 
sky. 



This paper is organized as follows. We discuss Stokes 
parameters in §2. Interferometric CMB polarization mea- 
surement on spherical sky is discussed in §3. In §4, we show 
visibilities are linearly weighted sum of spin ±2 spherical 
harmonic coefficients. In §5, we show how spin ±2 spherical 
harmonic coefficients can be determined from visibilities. In 
§6, computational feasibility is discussed. In §7, reconstruc- 
tion results from simulated observations are presented. In §8, 
the summary and discussion are given. In Appendix A, we 
discuss methods to facilitate computation of linear weight 
for o±2,im- In Appendix B, the reconstruction results with- 
out noise is presented. 
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2 ALL-SKY STOKES PARAMETERS 



3 INTERFEROMETRIC MEASUREMENT 



There are Stokes pa rameter Q and U, which describe th e 
state of polarization l|Kraudl 19861 ; iRohlfs and Wilsonll2003l) . 
Since Thomson scattering does not generate circular polar- 
ization in early Universe, circular polarization state V is 
not considered here. In this pap er, we follow the p olariza- 
tion convention of the healpix l|Gorski et al . 1 120051 ) . which 
differs from the definition of the International Astronomi- 
cal Union. In all -sky analysis, these are measured in refer- 
ence to (eg,efi) (|Zal darriaga 1998al; [Zaldarriaga and Seliakl 
1 19971 ). eg and are unit vectors of a spheric al coordinate 
system and given by (jArfken and Webe: 



sphcrica 



eg — i cos 9 cos (f> + j cos 6* sin — k sin#, 
= — i sin^+j cos</>. 

Stokes parameters Q and U are as follows: 

Q = (Eg — E 2 ^ , 
U = {2E e Ei) , 



(1) 

(2) 



where (. . .) indicates time average. Q and U transform un- 
der rotation of an angle %[) on the plane perpendicular to 
direction n as 



<3'(n) = Q(n)cos2^ + £/(n)sin2V>, 
U'(h) = -Q(n)sin2?/' + £/(n)cos2V>, 



(3) 
(4) 



with which the following quantities can be c onstructed 
l|Zaldarriaga and Selia5ll997l ; IZaldarriagalll998ah : 

(Q±if/)'(n) =e =F2iV '(Q±i*7)(n). (5) 

For all-sky analysis, Q an d U are expanded in terms of 
spin ±2 spherical harmonics (|Zaldarriaga and Seliakl 1 19971 ; 
IZaldarriagalll998al rbT) as follows: 



Q(n)+iC/(n) = ]T 2 Yim(n), 

/ ,tn 

Q(h)~\U(h) = ^ a- 2 ,i m -2Y lm (n), 



(6) 

m 



a 2 .im is related to a 2 .im by a_ 2 , im = (™l) m a2, i-m 
l|Zaldarriaga and Seliakl Il997l. Spin ±2 spherical harmon- 
ics have following forms: 



= \l^^-lF 1 ,im(0)+F 2Jm (e)]e im t (8) 



2 Y lm {9,< 



21 + 1 
4n 



[F Um (9) - F 2 j m {9)]e lm,p , (9) 



where Fijm and F^.i m can be computed in terms of Legendre 
functi ons as follows ( Ka mionkowski et alJI 19971 : IZaldarriagal 
Il998al ): 



/ (I — 2) 1(1 — m) ! . . cos 8 „™ . 



,1-m' . 1, 
sin 2 



<— + ^1(1 ~l))Pr (cos 9)], (10) 



(l + 2)\(l + m)\ sin 2 9 ' 
-(I- 1) cos 9 Pr (cost 



(11) 



The discussion in this section is for an ideal interferometer. 
An interferometer measures time-averaged correlation of two 
electric field from a pair of identical apertures positioned at 
ri and at r 2 . The separation, B = ri — r 2 , of two apertures 
is called th e 'baseline' and the measured correla tion is called 
'visibility ' (|Lawsonll2006l ; lThompson et al.ll200ll) . Depending 
on the instrumental configuration, visibilities are associated 
with (El - E 2 y ), (2E x E y ) and (El - E 2 y ±\2E x E y ) respec- 
tively, where x and y are axes of the polarizer frame. As dis- 
cussed in Sj2] Stokes parameters at angular coordinate (6,<f>) 
are defined in respect to two basis vectors eg and e^. Con- 
sider the polarization observation, whose antenna pointing 
is in the direction of angular coordinate (9a,4>a)- The po- 
larizers and baselines are assumed to be on the aperture 
plane. Then, the global frame coincides with the polarizer 
frame after Euler rotation (<j>A, 9a, VO on l ne global frame, 
where ip is the rotation around the axis in the direction of 
antenna pointing. Most of interferometer experiments for 
the CMB observation employ feedhorns for beam collection. 
After passing through the feedhorn system, an incoming off- 
axis ray becomes on-axis ray. Then, the basis vectors eg and 
e<p of the ray after the feedhorn system are related to the 
basis vectors e x and e y of the polarizer frame as follows: 



e^ + iey^e (&e A + ie0 A ) 
where $ is given by 



,i(*-V0/ a 



(eg + ie^ 



(12) 



4> 



tan 



+ tan 



sin 9 cos 9a cos(0 — 4>a) — cos 9 sin 9a 
sin 9a sin(0 — 0a) 



- sin 9 cos 9a + cos 9 sin 9a cos(0 — <f>A) 

Refer to Appendix [A] for the details on the derivation of <f> 
With Eq. 1121 we can easily show that 



{El - E 2 } + i<2£;, E y ) = e 



i(2i/>-2$) 



((E 2 e-El)+i{2EeE^)). 



With the employment of linear polarizers, the visibilities are 
associated with (El — E 2 ) or (2E x E y ), and are as follows: 

V Q , = f(v) J dfL4(n,n A ) (13) 

xRe [e- i(2 ^^ 2,s(fl)) ((3(n) + if/(n))] e i27ru ' fl , 
Vv> = f(v) J dfM(n,n A ) (14) 

xlm [e- iW - 2 * (fl,) (Q(n) + iU(h))} e i2nu A , 

where Aa is the unit vector in the direction of antenna point- 
ing and f(v) is the frequency spectrum of the CMB polar- 
ization. □ With the employment of circular polarizers, the 
visibilities are associated with (El — E y ± i 2E x E y ), and are 
as follows: 

V RL = f(y) J d04(n,riA), (15) 

X [Q(n) + i[ /( fi )] e H2 I ru.fl-2V. + 24»(n)) _ 

Vlr = f(v) I AQA(h,h A ), (16) 



8B(v,T) I 

ST \t=T 
To is the CMB monopole temperature 



1 m 



where B(u, T) is the Plank function and 
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x[Q(n) - i[ /( fi )]e i ( 2 ' ru ' fl + 2 *- 2 *( fl )) i 

where R and L stand for right/left circular polarizers. 

As in the following, Vqi and Vjji are linear combinations 
of Vrl and Vlr, and vice versa. 



-(Vrl + Vlr), 



V v > = -^(Val- Vlr), 
Vrl = Vq,+iVw, 
Vlr = Vq/ -iVc,/, 



(17) 

(18) 

(19) 
(20) 



4 VISIBILITY AS THE LINEAR SUM OF 
SPHERICAL HARMONIC COEFFICIENTS 

With Eq. [S] and [7] CMB visibilities can be expressed as a 
linearly weighted sum of a±2,i m in following ways: 



Vqi — — J^(a2,im&2,im + «-2,im &-2,im ) , 



(21) 



/ . m 



Vxjt = — — ^ ^ (Q2,Zm^2,Zm — 0>—2,lmb-2 7 lm)-> (22) 

(23) 
(24) 



Vrl = > y a 2,imi>2,lm, 
VLR = ^ ' a-2,lrqb-2,lm, 



where 

fe 2 ,im = y dvf{y) J dfL4(n, n A ) 

Xe i(2.u i .n- 2 ^ +2 *Cfi)) 2 y im(fi)i 

6-2,im = y duf(v) J dfL4(n, ua) 

xel (2xu, A+ 2*-2#(n))_ 2 y im(fi)i 



(25) 



(26) 



All the instrumental and configurational information arc 
contained in 6±2,im- As seen in Eq. 121 1 [2211231 and 1241 6±2, im 
are linear weights for a±2d m - As seen in Eq|25land[26l 6±2,im 
have distinct values, which depend on its spherical harmonic 
number I, m and the observational configuration such as an- 
tenna pointing and baseline. 



5 DETERMINATION OF INDIVIDUAL 

SPHERICAL HARMONIC COEFFICIENT 

An interferometer is sensitive to multipoles of a range lo ^ 
I ^ li. °~t 1 is given by 27ru, where u is a baselinelength 
divided by a wavelength. The width of the range, h — lo, 
depends on the window function, which is the square of the 
beam function in spherical harmonic space. When the inter- 
ferometer is sensitive to multipole of a range lo I ^ 1%, 
there are (h + l) 2 — 1% spin ±2 spherical harmonics in 
the range. Spin ±2 spherical harmonic coefficient a±2,im 
(lo ^ I ^ h) can be determined by fitting them for vis- 
ibilities of various antenna pointings and baseline orienta- 
tions. For coding convenience, we split visibilities, b±2.im, 
and a±2,im into real and imaginary parts. We enumerate 



real and imaginary parts of visibilities, b±2,im, and a±2,i m 
in matrix notation as follows: 

V = ba. (27) 
Likelihood function is given by 

£ = , w! ^ exphl(A-ba) T N- 1 (A-ba)], 
(27r)"(detN)2 2 

where n is the number of visibilities, A is a visibility data 
vector, N is a noise covariance matrix. The likelihood func- 
tion is maximum at 



a= (b T N" 1 b)" 1 b T N" 1 A. 



(28) 



Eq.[28]is reduced to a = b X A if b is square and b is invert- 
ible. The covariance of estimation error is 

(Aa Aa T ) 

= <(b T N- 1 b)- 1 b T N- 1 A JV {(b T N- 1 b)- 1 b T N- 1 A JV } T >, 
= (b T N _1 b) _1 , (29) 

where Ajv is the noise of a visibility data vector and 
(AjvAjv T ) = N. The E and B decomposition modes can 
be determined as follows: 

(30) 



= -(0 2 ,im + (-l) m d2,i -m)/2 
= -((-l) m al 2 ,i-m + a_2,Im)/2, 



CbB,lm ~ i(«2, im — ( — l) m ffl2,! -m)/2 

= i((-l)" l al 2 ,i-m - a-2,!m)/2. 
The variance of clei™. and a,B,lm estimation error are 

/ I ) ) / 

= T((Aa ±2 , !m Aa±2,im) + (Aa± 2 ,!-mAa± 2 ,i-m) 



+ (-l) m 2Re[{Aa ±2 , im Aa4 2 



*>])> 



(31) 



(32) 



(33) 



= -({Afl ±2 ,lmAa±2,Im) + (Aa ±2 ,j - m Ao± 2 ,l -m) 

-(-l) m 2Re[(Aa ±2 , im Aai 2 , ! - m )]), 

where the variance and covariance of Aa± 2i ; ± m are given by 
Eq.ESl 

6 SCALING OF COMPUTATIONAL LOAD 

As shown in previous section, a 2 ,!m (la ^ I ^ h) is deter- 
mined by 

a= (b T N" 1 b)" 1 b T N" 1 A. (34) 

a is the vector of length m, b is a n x m matrix and N is a 
n x n matrix, where n is the number of visibilities and m is 
the number of 02, i m , which is (h + l) 2 — /§. Unlikes the instru- 
mental noise of a single dish experiment, the noise covariance 
matrix fo r interferometric observations can be assumed to be 
diagonal llPark et al.ll2003h . Since inverting matrix is 0(Af 3 ) 
ijPress et al.lll992h while inverting diagonal matrix is 0(M) 
operation, Eq. [34] is a process of C^m 3 ). Computing b±2,im 
is small computational load, compared with computing Eq. 
1341 The method to compute b±2,i m fast is presented in Ap- 
pendix[B] Rough estimate by scaling our simulation in Sj7]to 
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higher multipole (I ~ 1500) says it will take roughly ~ 2500 
days by Pentium 4 (2Ghz) system. With ultra performance 
of computers such as SGI or IBM, determination of a,2,im 
over high multipoles by this formalism is computationally 
feasible. 



7 SIMULATED OBSERVATION 



We used the CAMB (|Lewis et all l2000h to compute the 
power spectra of ACDM and tensor-to-scalar ratio (r = 
0.3). a,E,tm and a,B,tm sets are drawn from the CAMB power 
spectra. With these aE,im and asdm, we have generated the 
simulated CMB Q and U maps by 

a±2,im ±2Yi m (n). 

l,m 

With the Q and U maps, Vrl and Vlr were simulated by 
numerically computing the following: 

Vrl = noise + / duf(u) I dflA(n, &a) 



x[Q(A) + iU(&)]e 



i(27TUi-ft-2i£-l-2*(ft)) 



Vlr = noise + / duf(u) I dflA(h, ua 



:[Q(n) 



We a ssumed the sensitivity of the Planck at 30GHz (|Tauberl 
2000): 13mJy (7.64 /iK) for noise, though nature of instru- 
ments are different. The observational frequency was as- 
sumed to be 30 - 31 GHz with 30° FWHM Gaussian pri- 
mary beam0. Total visibilities is ng A x n^ A x n^ u , where the 
number of fields is ng A x n$ A . For each field, baseline orien- 
tations are assumed to be n/n^^ k, where k — 1,2, ■ ■ ■ n^^. 
This can be achieved by building a feedhorn array of a n^ u 
fold rotational symmetry. The fields of survey are assumed 
to have angular coordinate (tt i/ng A ,2n j/n$ A ), where 
i — 1,2, ■ ■ ■ ng A and j — 1, 2, ■ ■ ■ n^ A . As seen in Eq. IBll and 
IB21 b±2,i m depends on antenna pointing (0a,4>a) and base- 
line orientation <f> u . In simulated observation, we assumed 
ne A = n 4>A = = n 1/3 > ((h + if - Zg) 1/3 for the 

variation of (9a, (f>A,<fiu)- 

Baselines of length B\ = 6 [cm], B2 — 9 [cm], B3 = 
12 [cm] and B 4 = 15 [cm] with 30° FWHM beams were 
assumed. The corresponding window functions are shown 
in Fig. [1] A window function corresponding to the longest 
baseline is shown at rightmost. The interferometers of Bi 
are sensitive to the multipole range 29 ^ I ^ 48, those of 
B 2 to 48 < I < 67, those of B 3 to 67 < I < 86 and those 
of B4 to 86 ^ I ^ 106. We chose the multipole range lo 
and h to be the first multipoles where the window function 
drops below 1% of its peak value. With such cutoff, there 
exists error from residuals, which contributes to estimation 
error. We chose ng A — n^ A — n^> u = 20, which makes the 
total number of visibilities (Vrl and Vlr) about four times 
the number of the spherical harmonics to be determined. So 
the number of constraints is about four times the number of 



2 The conclusion of this paper is not affected by the shape of the 
beam function as far as the window function corresponding to the 
beam function is not non-zero over infinite number of multipoles. 
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Figure 1. Window Function (normalized to its peak.) 




Figure 2. estimated Rc[ag ; m ] (m=29) 



unknowns, which is necessary in the presence of noise and 
residual error. 

With Eq. 1281 we have estimated spherical harmonic co- 
efficients a,2,im (29 ^ I ^ 48) from Bi visibilities, (48 ^ I ^ 
67) from B 2 , (67 < I < 86) from B 3 , and (86 < I < 106) 
from B4. From the estimated a2,; m , we have obtained CLE,im 
and as.im via Eq. 1301 and |3T1 Estimated aE,lm and an,im 
( m — 29) are shown together with the input a,E,im and 
o-B,im value in Fig. [2] [3] [4] and [5] These are representative 
of other m modes. The 1 — a errors via Eq. [29] are indicated 
by vertical error bars in Fig. [21 El SI and[S] As shown in Fig. 
[1] the window function of the simulated observation have 
several troughs, where 1 — a errors are expected to be large. 
It is seen that 1 — a errors are significant on the multipoles 
corresponding to the troughs of the window function. By 
averaging the magnitude of 1-a error on a,E,im and a,B,im, 
we found that they are in same magnitude within 1%. It is 
not suprising, considering the expression for the variance of 
AaB,i m and Aa,B,im, which are shown in Eq. 1321 and 1331 
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Figure 3. estimated lva[a^ ; m ] (m=29) 
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Figure 5. estimated lm[a B tm ] (m=29) 



8 DISCUSSION 

Visibilities associated with the CMB polarization can be ex- 
pressed as a linear sum of spherical harmonic coefficients 
0±2, ira of the CMB polarization. The linear weight for a±2,i m 
depends on the observational configuration, and spherical 
harmonic number I, m. Since an interferometer is sensitive 
over a finite range of multipoles. The spherical harmonic co- 
efficients (lo ^ I ^ h) can be determined by fitting a±2,i m 
for visibilities of various observational configuration. Once 
a±2,im are determined, a,E,im and a,B,im are easily obtained 
via Eq. l30l and l31l The linear weights 6±2,im, which map vis- 
ibilities to spherical harmonic space, can be computed fast 
with the aid of the methods discussed in Appendix [B] The 
best-fit value of a2,im for given visibilities can be found via 
Eq.[28] It is C(m 3 ) process, where m is the number of a,2,i m 
to be determined. Scaling the time taken for the simulated 
observation says this formalism is computationally feasible 
for interferometric observation up to multipoles as high as 
I ~ 1500. 

Since the formalism introduced in this paper is devel- 
oped for a satelite-based interferom etric observations o f all- 
sky polarization such as the EPIC (|Timbie et al.|[200r3 ). the 
antenna pointings in the simulated observation are made 
over full-sky. Even when antenna pointings are made within 
a fraction of sky, the matrix b in Eq. [23 does not become 
singular, as far as the fraction of sky is big enough in relative 
to the angular scales the interferometer is sensitive to. 

a±2,im determined by Eq. [25] contain foreground con- 
tamination like pixellized maps. With different frequency 
spectral behavior of foregrounds from the CMB, the fore- 
ground contribution can be separated with the multi- 
frequency data from the CMB down to residual level, which 
is limited by frequency spec tral incoherence and knowledge 
of polarized foregro unds l|Tegmark and Efstathioul Il996l ; 
iTegmark et~aT]|2000l ). 

There are several sources for E/B mode mixing. Alias- 
ing due to finite pixel size leads to E/B mode mixing at high 
multipoles and limited sky coverage does at low multipoles. 
Interferometer observations enable targeting a specific range 
of multipoles. E/B mode mixing at low multipoles due to 
limited sky coverage can be made insignificant by designing 
interferometers to be insensitive to anisotropy at low mul- 
tipoles. Since the formalism reconstructs a,E,im and aB,im 
directly from spherical harmonic space, E/B mixing due to 
pixellization and ambiguity of E/B mode over the mosaic 
are insignificant. 

We choose the multipole range lo and h to be the first 
multipoles, where the window function drops below 1% of 
its peak value. There are residual contribution from a±2,im 
(I < lo, I > h). These residuals are another source of estima- 
tion error in addition to instrument noise. We can modify 
noise covariance matrix of Eq. [28] to include y], C[W|.„, 
where Ci and W^y is power spectra and window func- 
tions at out-of-bound multipoles. It improves the estima- 
tion error due to residuals by giving more weights to the 
visibilities of less contribution from residuals. But it in- 
creases computational load, by making a total noise co- 
variance matrix non-diagonal, while an instrument noise 
covariance matrix is diagonal t o a good approximation 
l|M.P.Hobson and Maisingerll2002l ). 

We have presented a formalism to reconstruct spherical 
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harmonics of the CMB polarization directly from interfer- 
ometer observations. The formalism takes advantage of the 
fact that an interferometer directly probes the Fourier com- 
ponents of sky pattern, and the relation between a Fourier 
component and spin ±2 spherical harmonics. 
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APPENDIX A: CMB POLARIZATION BASIS 
VECTORS AND ANTENNA COORDINATE 

In all-sky analysis, the CMB polarization at the angular 
coordinate (9, <f>) are measured in the local reference frame 
whose axises are (eg, e,/,, e r ). Let's call this coordinate frame 
'the local CMBP frame' from now on. Consider the polar- 
ization observation of antenna pointing (9a,4>a). A global 
coordinate frame coincides with the antenna coordinate by 
Euler rotations R !/ (6a) H. z (4>a). Since a global coordinate 
frame coincides with the local CMBP frame by Euler rota- 
tions Ry(9) R z (</>), the Euler Rotations R z ( 7 ) R M (/3) R z (a) 
coincides the antenna coordinate frame with the local 
CMBP frame, where R z ( 7 ) R H (/3) R z (a) R. y (9 A ) R z (>a) = 
Ky(9)K z (<j)). Therefore, the local CMBP frame is in ro- 
tation from the antenna coordinate by the Euler angles 
(q, (3, 7) as follows: 



a = tan 



sin 9 sin(^> — 4>a) 



sin 9 cos 9a cos(0 — </>a) — cos 9 sin 9a 
(3 = cos -1 [cos 9 cos 9a + sin 9 sin 9a cos((j> — cj>A)] ■ 



tan 



sin 9a sinf 



- sin 9 cos 9a + cos 9 sin 9a cos(0 — 4>a) 



where the Euler angles (a, f3, 7) can be obtained from 



R z ( 7 )R„(/3)R z (a) = R H (0)R z (»R; 



iR-'^a). In 



most CMB polarization experiments, where polarizers are 
attached to the other side of feedhorns, incoming rays go 
through polarizers after feedhorns. After passing through a 
feedhorn, an incoming off-axis ray becomes an on-axis ray. 
Then the local CMBP frame of the ray after the feedhorn 
system is simply in azimuthal rotation a + 7 from the an- 
tenna coordinate. Therefore, $ in Eq. 1 131 is 



$ = tan 



+ tan 



sin 9 sin(0 — 4>a) 



sin o cos Ha cos( 



cos # sin Ha 



sin 9a sin((p — <f>A) 



+ cos 9 sin 9a cos((p — 4>a) 



APPENDIX B: COMPUTING LINEAR 
WEIGHTS 

b±2,im needs be computed to determine a% t i m {lo ^ I ^ h) 
via Eq. 1281 It can be computed in the baseline coordinate 
where a baseline coincides with the x axis of the coordinate. 
In computing b±2,i, n in the baseline coordinate, spin ±2 
spherical harmonics, which are defined in the global coordi- 
nate, are related to spin ±2 spherical h armonics in the base- 
line c oordinate with a rotation matrix ( Cha llinor and Lewis! 
[20051 ). Let's choose Galactic coordinate as the global ref- 
erence coordinate for the CMB. Consider the polarization 
observation, whose antenna pointing is in the direction of 
Galactic coordinate (£,b). The baseline is assumed to be on 
the aperture plane. Then, the baseline coordinate is the co- 
ordinate system rotated from Galactic coordinate by Euler 
rotation 1Z(£, \ — b, U ), where ^ u is the rotation around the 
axis of antenna pointing. Computed in the baseline coordi- 
nate, b±2,im are 



02,lm — 



duf(y) I d(9') sin 9' A(9') 
Jo 

2tt 



i>2n 

J d4>' ^D^ m (£,|-6,^ u ) 2 Y lm ,(e' ',</>') 



xe 

— e 



i(27rusine' cos </>' -2(i/>-0 u )+2</>' ) 



- i2 V^ + i) 



x / d(cos9')A(9')[F Um ,(9') + F 2 , lm ,(9')] 



-2,1m 



x / dv f(u) J m / +2 (27r sin 9'u), 



dvj(y) \ d(cos(9')A(60 



(Bl) 



xe 

— e 



^ m' 
i(27rusin# cos <f> - 

i2i '^n(2l + l) 



b,4>v) -2Y lm >{9' ,4>'\ 

\(2ttu sin 9' cos <p' +2(i> - 4> m ) — 2ip' ) 
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x / d(cos 0^(00 [Fi^O-i^w^')] 



x / du f(u)J m i _ 2 (2tt sin 9' u), 



(B2) 5- 



where J n is the nth order ordinary Bessel function. It turns 
out that computing b±2. im via Eq. IB1I and IB2I suffers from 
serious numerical precision problem especially for high mul- 
tipole (I > 60), which are due to machine floating-point 
rounding error occurring in the multiplication with the rota- 
tion matrix D l m , m (1Z). Besides the numerical precision prob- 
lem, the huge time required for computing the rotation ma- 
trix D l m , m (lZ) makes Eg. IBll and Eg. IB2l lose most of merits 
gained by the availability of analytic integration over az- 
imuthal angle. For these reasons, in the simulated observa- 
tion of 33 we computed 6±2, im in a fixed CMB frame with 
Eq. [25] and [26] In Eq. [25] and 1261 we have rearranged the 
order of integration and replaced the integration over con- 
tinuum with sum over finite elements, which are as follows: 



62,1 



21 + 1 -\2ih 

— e 

4% 



A9A<j>Av 



(B3) 



^sm0i(Fi,i m (0i)+.Fvm(0i)) 



x^e 



m<j>j i2*(n) ' 



k 



f(v k )A(h, n A )e 



i(2iru fc -n) 



ll±±e i2 ^A9A^Av 

+K 



(B4) 



x^sin0 < (F 1>lm (0 j )-.F2,j m (0 < )) 
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\ A imiA. — i2<&(ri) \ A £ r \ a f *• - \ i(27rui.-n) 
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Figure CI. estimated Re[a^ i m ] (m=29) in absence of noise 



form (|Press et all Il992l ; iMuciaccia et al.lH997T ). Since it is 
easiest to carry out Fast Fourier Transform on data of 
number which is a power of two, A<f> was chosen to be 
360°/(2 riog2(2!l)1 ), where |~] denotes the smallest integer 
larger than or equal to the argument. Choosing the optimal 
size for A0 and and using Fast Fourier Transform enables 
the numerical computation of b±2,i m feasible in a reasonable 
amount of time even for interferometers of high u. We set 
the value for the integration cell to be zero and skips com- 
puting the rest of terms, when the beam function A(n, ua) 
for the integration cell is smaller than .1% of its peak value. 
In integration over bandwidth, f(yk)A(n, nA)e l2,rB ~^" n , 
an interference term of index k is computed from a term of 
index k — 1 as follows so that we can avoid computing time- 
consuming trigonometric function for each index k: 



Legendre functions in i 7 i l i m (0) and F2,i m (9) are com- 
puted with the fo llowing recurrence relation for x = cos 8 
ijPress et al.lll992T ): 

(l-m)Pr = x(2l-l)P^ 1 -(l + m-l)P{1 2 , 
PZ = (-l) m (2m-l)!!(l- a; 2 ) m - 2 , 
P™ +1 = x{2m + l)P™. 

When an interferometer is sensitive to multipole range lo ^ 
I ^ h, b±2,im of / up to h should be computed. As shown 
in Eq. [8] Eq. [9] Eq. [10] and Eq. [11] ±2 Y tm is the product 
of Legendre functions and e !m< ^. Legendre functions of mul- 
tipole I varies on angular scale dow n to 9 « 180° /I. With 
Nyquist-Shannon Sampling theorem l|Press et al.ll 19921 ). the 
integration over 9 should be done with A9 < Since 
tgim^j cons j s ting of ±2Yi rn is periodic over 4> = 360° /m, the 
integration over cj> should be also done with A(f> < ^p-. 

The summation of index j below, which is part of Eq. 
IB3l and lB4l is equivalent to discrete Fourier Transform: 



e im ^ 



. Discrete Fourier Transform, which is the process of 0(JV 2 ), 
can be carried out in 0(A/"log 2 A/") with Fast Fourier Trans- 



2tt Av . . 2tt Av i27r ^i B .fi 

(cos( B-n) + ism( B ■ n))e <= , 

c c 

where the computed value of cos( 27r e Al/ B-n) and sin( 2,r Al/ B- 
n) are repeatedly used. 



APPENDIX C: ESTIMATION IN THE 
ABSENCE OF NOISE 

Estimated a 2 ,i m (29 ^ Z ^ 106, m = 29, 0, -29) in ab- 
sence of noise are shown together with the input a2,; m value 
in Fig. [CI] \C2\ [C3] and We assumed the same con- 
figuration with the simulated observation in !j7] except for 
the absence of noise. Small discrepancies between estima- 
tion and the input values, in spite of no noise, are attributed 
to residual error. As mentioned in ij7] the residual error re- 
sults from the contribution of 02,i m in the multipoles outside 
the cutoff region, since we determined a2,; m only over the 
multipoles where the window function is greater than 1% of 
its peak value. Some features of methods to facilitate 62, im 
computation, which are discussed in Appendix [B] sacrifice 
the numerical precision, which also contributes to the dis- 
crepancies. 
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Figure C2. estimated Im[<2E ; m ] (m=29) in absence of noise 




Figure C3. estimated Re[aa ; m ] (m=29) in absence of noise 
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Figure C4. estimated lm[a B im ] (m=29) in absence of noise 



